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ABSTRACT 

The present paper is the last of a series studying the first-order Fermi acceleration processes at 
relativistic shock waves with the method of Monte Carlo simulations applied to shocks propagating 
in reali stically modeled turbulent magnet ic fields. The model of the background magnetic field struc- 
ture of iNiemiec fc Ostrowskil l)20f)4 1200(11) has been augmented here by a large-amplitude short-wave 
do wnstream component, i mitati ng that generated by plasma instabilities at the shock front. Follow- 
ing [Niemiec^^Jpstrowski (2006), we have considered ultrarelativistic shocks with the mean magnetic 
field oriented both oblique and parallel to the shock normal. For both cases simulations have been 
performed for different choices of magnetic field perturbations, represented by various wave power 
spectra within a wide wavevector range. The results show that the introduction of the short-wave 
component downstream of the shock is not sufficient to produce power-law particle spectra with the 
"universal" spectral index 4.2. On the contrary, concave spectra with cutoffs are preferentially formed, 
the curvature and cutoff energy being dependent on the properties of turbulence. Our results suggest 
that the electromagnetic emission observed from astrophysical sites with relativistic jets, e.g. active 
galactic nuclei and gamma-ray bursts, is likely generated by particles accelerated in processes other 
than the widely invoked first-order Fermi mechanism. 

Subject headings: acceleration of particles, cosmic rays, gamma rays: bursts, methods: numerical, 
relativity, shock waves 



1. INTRODUCTION 

Realistic modeling of first-order Fermi cosmic-ray ac- 
celeration at relativistic shock waves is a difficult task 
due to the strong dependence of the resulting par- 
ticle spectra on the essentially unknown local con- 
ditions at such shocks. Therefore, progress in this 
field mostly arises from the increasingly more realistic 
MHD conditions that are assumed to apply near the 
shocks. For mildly relativistic shocks in the test par- 
ticle approach the advances that have been made in 
the past twenty years involve a semianalytic solution 
of the particle pitch-angle diffu sion equation at parallel 
shocks ijKirk fc Schneider] Il9 87). the effects of different 
shock compressions in parallel shocks ({Heavens fc Drurvl 
1988), the impact of an obliq ue mean magnet i c fiel d 
for subluminal configurations (Kir k fc Heavens) 11989^1 . 
and the effects of oblique field configurations at 
sub- and superluminal shocks considered with an in- 
creasing degree of complication, appl ying successive! 
more realistic physical approxima t ions (|O strowsk 3ll99 
19931 IBednarz fc Ostrowskil Il99fit INiemiec fc Ostrowsk 
2004). For ultrarelativistic shocks the notion of 
a "universal" spectral inde x was first described by 
IBednarz fc Ostrowskil l|1998fl and later found in a va- 
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riety of studies llCallant fc Ac hterberdlT999t iKirk et all 
1 20001 lAchterberg et alJl2001t iLemoine fc Pelletierl 12003 

Ellison & Double 200 4), b ut more recent studies 
(jNiemiec fc Ostrowskil l2006fl indicate basic problems 
with first-order Fermi acceleration, in particular in gen- 
erating wide-energy power-law spectra. The latter work 
also contains a more systematic discussion of the subject 
literature. 

On the other hand, there have been attempts to self- 
consistently derive the electromagnetic shock structure. 
One very simplified approach applied modeling of non- 
linear effects in the first-order Fermi acceleration at the 
shock (e.g., Ellison & Double 2002). The limitations of 
this technique arise from the approximate treatment of 
particle scattering near the shock and its simple scaling 
with particle energy instead of a realistic accounting for 
the microphysics of the process. A more realistic mi- 
croscopic and self-consistent description of the collision- 
less shock transition ma y be afforded by Particle-In-Cell 
(PIC ) simulations (e.g. . iHoshino et al. [1992]: [Sily^et^l " 
2003; Nishikawaet all 120031 l2005at Frederiksen et al 
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which follow the formation of a relativistic collisionless 
shock starting with processes acting on the plasma scale. 
In these simulations one typically observes that the gen- 
eration of magnetic and electric fields is accompanied by 
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an evolution of the particle distribution function from the 
background plasma conditions up to highly superthermal 
energy scales. Although providing a wealth of informa- 
tion on the collisionless shock structure, the PIC simu- 
lations are still very much limited in dynamical range, 
in particular for cosmic-ray particles whose energies are 
many orders of magnitude larger than the plasma parti- 
cle energies. 

The present paper i s a continua tion of the studies of 
iNiemiec fc Ostrowskil l|200l H)06). Here we attempt to 
incorporate the results of the PIC simulations concern- 
ing the generation of magnetic fields at relativistic shocks 
into the test particle Monte Carlo simulations of cosmic- 
ray particle acceleration. As in the previous papers, 
we study ultrarelativistic shock waves propagating in a 
medium with a magnetic held that is perturbed over a 
wide range of macroscopic scales. However, we also con- 
sider an additional short-wave isotropic turbulent held 
component downstream of the shock, analogous to the 
shock generated turbulent holds revealed by the PIC sim- 
ulations. We treat the amplitude of the short-wave com- 
ponent as a model parameter and study the effects of 
the small-scale perturbations on the particle spectra and 
angular distributions derived in the previous work. 

In what follows, c = 1 is the speed of light. The inte- 
gration of particle trajectories is performed in the respec- 
tive local plasma (upstream or downstream) rest frame, 
and we use index "1" ("2") to label quantities provided 
in the upstream (downstream) frame. We consider ul- 
trarelativistic particles with p = E. We use dimension- 
less variables, so a particle of unit energy moving in an 
uniform mean upstream magnetic field B$ has the unit 
maximum (for p± = E) gyroradius r g (E = 1) = 1 and 
the respective resonance wavevector is k res (E = 1) — 2tt. 

2. DETAILS OF THE MODELING 

2.1. Magnetic Field Structure 

As in INiemiec fc Ostrowskil i|2004 HflOfi) we assume a 
relativistic shock wave to propagate with velocity i*i (or 
the respective Lorentz factor 71) in the upstream medium 
with a uniform magnetic field Bo = -Bo,ii inclined at an 
angle ipi to the shock normal (along Ui direction), su- 
perimposed on which are isotropic 3D magnetic field per- 
turbations. These perturbations are characterized with a 
power-law wave power spectrum F(k) ~ kr q that is de- 
fined in a wide wavevector range {k\, m im ki^ max ), where 
ki,min = 0.0001 and k\^ max = 10. Specifically we con- 
sider the wave spectral indices q = 1, describing the flat 
power spectrum, and q — 5/3 for the Kolmogorov dis- 
tribution. The integral power of these perturbations is 
given by the (upstream) amplitude 

SB = ^ J 1,m °* F(k)dk, (1) 

and SB /Bo is one of our model paramete r s. As described 
in detail in INiemiec fc Ostrowskil 1)2004 I2006JL the up- 
stream field perturbations are modeled as the superpo- 
sition of static sinusoidal waves of finite amplitude. The 
downstream magnetic field structure, including the tur- 
bulent component, is derived as the shock-compressed 
upstream held, and hence the downstream turbulent 
magnetic field is naturally anisotropic. Our method of 
applying the hydrodynamical shock jump conditions (for 



the electron-proton plasma. iHeavens fc Drurvlf l988) pre- 
serves the continuity of the magnetic field lines across the 
shock. In what follows, we refer to the field component 
described above as the large-scale (or long-wave) back- 
ground field. 

The other physical magnetic held component consid- 
ered in the present study is the short-wave turbulence, 
assumed to be a result of kinetic magnetic field gener- 
ation processes acting at the shock. This component, 
taken for simplicity to be isotropic 1 and static, is imposed 
upon the nonuniform background magnetic held down- 
stream of the shock. The short-wave nonlinear field per- 
turbations are introduced with a flat spectral distribution 
in the wavevector range (10&2, max, 100&2, max), where 
the shortest downstream waves are the shock-compressed 
shortest upstream waves: 

k*2,^max — ki^max B . (2) 
72 

Here R = u\fui is the compression ratio in the shock rest 
frame, and 71 (72) is the upstream (downstream) shock 
Lorentz factor. The wavevector range of the small-scale 
turbulent component is chosen to be one decade in k 
apart from k2, ma x to separate the influence of these per- 
turbations on the low-energy particle motion from that 
exerted by the short-wave component of the large-scale 
background field. This choice also facilitates the use of 
a hybrid method for the calculation of particle trajec- 
tories (see §2.2). Note also, that the definition of the 
short-wave component depends on the Lorentz factor of 
the shock. Below, we refer to the short-wave magnetic 
field component also as the shock- generated turbulence 
( "sh" ) , to distinguish it from the large-scale background 
field, that exists in the upstream region and is only com- 
pressed upon passage through the shock. 

In our simulations we exclusively study the first-order 
Fermi acceleration process and neglect second-order pro- 
cesses. Therefore, the turbulent magnetic field compo- 
nents (both short-wave and large-scale background per- 
turbations) can be considered to be static in the respec- 
tive plasma rest frames, both upstream and downstream 
of the shock, and electric fields that may exist in the 
shock transition layer are neglected. 

2.2. Monte Carlo Simulations 

The implementation of short-wave turbulence into the 
Monte Carlo simulations forces us to dispense with 
the direct integration of the particle equations of mo- 
tion in the analytically m odeled magnetic field (see 
INiemiec fc Ostrowskil [2006) . Instead, we resume the 
derivation of parti cle trajectories with the hyb rid ap- 
proach proposed in INiemiec" fc Ostrowskil l)2004j) . Thus, 
particle trajectories are directly calculated from the 
equations of motion in the large-scale background mag- 
netic field only, whereas the trajectory perturbations due 

L Theoretical considerations byJMe^edev&LoetJ_^)99^_jMid 
numerical PIC simulations (e .g.. lSiL^^T^ni20nSTNTsru^awaet alJ 

2003; Frcdcrikscn ct al. 2004) show that the relativistic two-stream 
(Weibel-like) instability at the relativistic collisionless shock front 
leads to the generation of a strong, small-scale turbulent field down- 
stream, that is predominantly transversal (2D) and lies in the plane 
of the shock. Because the magnetic field in the nonlinear regime 
of the two-stream instability is sustained by the structure of the 
ion current channels, instabilities in the ion filaments (e.g., kink 
and/or firehose instability) should lead ultimately to the 3D tur- 
bulent magnetic field. 
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to the shock-generated small-scale turbulence are ac- 
counted for through a small-amplitude pitch-angle scat- 
tering term. Auxiliary simulations have been performed 
to determine the scattering amplitude distributions for 
various particle energies. In the scattering procedure, af- 
ter each time step At the particle momentum direction 
is perturbed by a small angle Af2. The time step itself 
scales with particle energy (or gyroradius), At oc E, and 
is chosen so that the condition cAt ^> X s h is always ful- 
filled (X s h is the wavelength of shock-generated perturba- 
tions), which means that the auxiliary simulations follow 
the particle scattering on the shock-generated turbulence 
well into the diffusive regime. Therefore, the scatter- 
ing amplitude distributions also scale linearly with the 
short-wave turbulence amplitude, and the mean scatter- 
ing angle variations are related to the particle energy as 
AflotE- 1 / 2 . 

An individual simulation run is performed as follows. 
We start with injecting monoenergetic particles (with 
the initial energy Eq — 0.1) at random positions along 
the shock front, with their momenta isotropically dis- 
tributed within a cone around the shock normal pointing 
upstream of the shock. In an initial simulation cycle, the 
injection process is continued until the required number 
of particles, N, has been selected — those that after being 
injected upstream and then transmitted downstream of 
the shock, succeed in recrossing the shock front again. 
Then the calculation of individual particle trajectories 
proceeds through all subsequent upstream-downstream 
cycles. In each cycle, a fraction of the particles escapes 
through a free-escape boundary introduced "far down- 
stream" of the shock, i.e. at a location from which there 
is only a negligible chance that particles crossing the 
boundary would return back to the shock. To the parti- 
cles remaining in the sim ulations the trajectory-split ting 
procedure is applied fsee lNiemiec fc O strowskil [2004D . so 
the number of particles remains constant in the acceler- 
ation process, but the statistical weights of the particles 
are appropriately reduced. The final spectra and angu- 
lar distributions of accelerated particles, derived in the 
shock normal rest frame for particles crossing the shock 
front, are averaged over many statistically different sim- 
ulation runs. 

3. RESULTS 

Because of the limited capabilities of present-day com- 
puters, the long-time nonlinear development of plasma 
instabilities leading to the generation of the short-wave 
turbulent downstream magnetic field component cannot 
yet be fully investigated, in particular for electron-ion 
plasma collision fronts. Therefore, the ultimate structure 
of the small-scale field and the global effectiveness of the 
generation mechanism remain uncertain. To estimate the 
role of the shock-generated turbulence in cosmic-ray ac- 
celeration at ultrarelativistic shocks, we introduce the 
small-scale wave component as isotropic 3D magnetic 
field perturbations and treat the amplitude of these per- 
turbations as a model parameter. 

The spectra of accelerated particles for obl ique super- 
luminal shocks are presented in Figures 1 1141 and those 
for parallel shocks are displayed in Figures El and |H1 
The amplitude of the short-wave component, 8B s h/ (B2), 
as provided in the figures, is measured in units of the 
average downstream perturbed magnetic field strength 



(B 2 ) = ({B 0y2 + 6B 2 ) 2 ) 1/2 . The spectra presented with 
solid lines has been derived in the turbulence model 
without the 8B s h term (m ost of them are presented in 
iNiemiec fc Ostr owski 2006). Some of the energetic par- 
ticle distributions follow a power-law in a certain energy 
range. In these cases linear fits to the power-law portions 
of the spectra are presented, and values of the phase- 
space spectral indices, a, are given (the equivalent num- 
ber spectral index a — a — 2). The spectral indices may 
help to make a quantitative comparison of the spectra, 
but care must be exercised, because the spectral indices 
depend on the energy range chosen for the fit on account 
of the curvature in most of the spectra. In the case of 
supcrluminal shocks, we use the same low-energy limit 
for the energy range selected for the fits. 

3.1. Oblique Superluminal Shocks 

Accelerated particle spectra for obliq ue superlumi- 
nal shocks are presented in Figures 11131 for if>i = 45° 
(ui/cos^i ~ l-4c), and in Figure E| for ipi = 90° (per- 
pendicular shock). One can see that increasing the am- 
plitude of the shock-generated turbulence leads to a more 
efficient acceleration with particle spectral tails extend- 
ing to higher energies. However, in all cases, in which 
fiB s h/{B2) ^ 1, the energetic spectral tails are convex, 
so the spectral index increases with particle energy In 
addition, all the spectra have cutoffs at an energy for 
which the resonance condition for interactions with the 
long-wave turbulence is fulfilled. These features result 
from the fact that the influence on particle trajectories 
of the shock-generated small-scale turbulence decreases 
with increasing particle energy. In our numerical ap- 
proach this corresponds to a reduction of the scattering 
amplitude AQ,(E) (see §2.2). 

Most published treatments of first-order Fermi accel- 
eration at ultrarelativistic shocks, that apply the pitch- 
angle diffusion approximati o n l|B e d narz fc OstrowskH 
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Ellison fc Double 
2005|), assume scattering condi- 
3 articl e energy. As 



tions that do not change with particl e 
discussed bv lOstrowski fc Bednard <f2002) . these authors 
show that if the small-angle scattering dominates over 
the possible influence of the oblique mean magnetic 
field component and/or the long- wave perturbations in 
shaping the particle trajectories, then power-law particle 
distributions with "universal" spectral index a u ~ 4.2 
may be formed. With the more realistic magnetic field 
model considered in the present paper, it is not possible 
to reproduce these results. However, particle spectra can 
approach a power-law form in a limited energy range 
near the particle injection energy, as demonstrated 
by the fits to the low-energy parts of the spectra in 
Figures [1141 The energy ranges for the fits are arbitrary 
selected but have the same low-energy bound. Note that 
in the figures the spectral indices, though depending 
on the background conditions, are always larger than 
a u . If one chose higher energies to derive a power-law 
fit to the spectra, one would obtain a higher spectral 
index on account of the convex shapes of the spectra. 
Convex spectra cou ld have been expected bas ed on 
the earlier studies of Bcdna rz fc Ostrowskil ((1998), who 
showed how an increased pitch-angle scattering angle in 
superluminal shocks leads to a flattening of the particle 
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Fig. 1. — Accelerated particle spectra at oblique superluminal 
shocks with 71 = 10, = 45°, and a flat wave-power spectrum of 
large-scale background magnetic field turbulence. The amplitude 
of the upstream field perturbations is SB/Bo = 1.0. The ampli- 
tudes of the downstream shock-generated turbulence, SBg^/ (B2), 
are given near the respective spectra. Linear fits to some spectra 
are also presented with dashed lines. The energy ranges selected 
for the fits have the same low-energy bounds to facilitate com- 
parison of the spectral indices a, given in italic. The spectrum 
plotted with a solid line applies to the case without short-wave 
shock-generated perturbations. The spectrum shown with filled 
dots is calculated with a fixed scattering term for E > 16, and the 
linear fit to the part of this spectrum above E > 16 is presented 
with dotted line. Some spectra are vertically shifted for clarity. 
Particles in the energy range (2-ir/kmax, 27r/fc m ; rl ) can satisfy the 
resonance condition k r es = 2ir/r g (E) for some of the waves in the 
background turbulence spectrum. 



spectrum to the limiting " universal" power law (see th e 
discussion of this effect in lOstrowski fc Bednard 120021. 
To further illustrate this behavior, we have recalculated 
the spectrum for 6B s h/(B2) = 80 in Figure ^ {filled 
dots) for a fixed scattering amplitude Afl(E) at particle 
energies E > 16 (AQ(E > 16) = const. < 1). Then 
the pitch-angle scattering term dominates downstream 
of the shock at high particle energies, which leads to a 
power-law spectrum formation over a wide energy range 
without a steepening or a cutoff, but with a > a u . 

In the case of superluminal shocks, the particle spec- 
tra calculated for a high amplitude of shock-generated 
turbulence, 8B s h/ (B 2 ) = 50 in Figures 111 II depend 
only weakly on the background large-scale magnetic field 
structure, in contrast to t he 5B s h = limit (see also 
iNiemiec fc Ostrows ki 2006). The acceleration process is 
only slightly more efficient in the case of a large ampli- 
tude of long-wave perturbations (5B/B n = 1.0 in Figs. ^ 
and [21; in which the spectral tails extend to marginally 
higher energies and show softer cutoffs, when compared 
to the cases with SB/B Q = 0.3 (Figs. Hand©. There 
are also only minor differences between the spectra de- 
rived for a flat spectrum of background perturbations, 
for which the wave power is uniformly distributed per 
logarithmic wavevector range, and for the Kolmogorov 
distribution, for which most power is carried by long 
waves. In the latter case, one finds power-law particle 
spectra in a slightly wider energy range than in the case 
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Fig. 2. — Particle spectra derived for oblique superluminal shocks 
for different amplitudes of shock-generated turbulence. For all 
spectra 71 = 10, = 45°, and a Kolmogorov power spectrum 
of background long- wave turbulence with SB /Bo = 1.0 has been 
assumed. The spectrum sho wn with the solid line is deriv ed in the 
limit SB 3 h = (see Fig. 2 in Niemiec & Ostrowski 2006). 
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Fig. 3. — Particle spectra at superluminal (71 = 10, Vi = 45°) 
shocks with SB/ Bo = 0.3. Filled and open symbols refer to results 
for a flat and a Kolmogorov wave-power spectrum of large-scale 
background field perturbations, respectively. The linear fit to the 
spectrum for the flat wave power spectrum is derived in an energy 
range with the same low-energy bound as used for the fits shown 
in Figs. [I]and|2] 



of flat distributions (compare Figs. and [21 and spectra 
in Fig. 0}, but the high-energy tails and the cutoff shape 
remain nearly the same for both types of the wave power 
spectra. 

The structure of the background large-scale magnetic 
field can, however, influence particle acceleration for 
smaller amplitudes of the short-wave turbulence. As 
one can see in Figure [21 the power-law part of the spec- 
trum for 8B s } l /(B2) = 20 and Kolmogorov- type long- 
wave turbulence spectrum, although steeper, continues 
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Fig. 4. — Particle spectra at perpendicular (71 = 10, jpi = 90°) 
shocks for the same parameter combinations as in Fig. y] 
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Fig. 5. — Normalized particle angular distributions as measured 
in the shock rest frame for superluminal shock waves with 71 = 10, 
ipl = 45°, and a Kolmogorov wave power spectrum of large-scale 
background magnetic field turbulence. Particles with cos 9 < 
are directed upstream of the shock. The distributions are calcu- 
lated by summing the quantity Mi/du^l + 0.005) in the respective 
cos 0-bins at every particle shock crossing, where 8 is the angle 
between the particle momentum and the shock normal, v x is the 
normal component of the particle velocity, and w is the statistical 
weight of the particle. Only particles with E > 1 are included, 
and low-energy particles with large weights contribute most to the 
distributions. The angular distribution presented with filled dots 
is derived without short-wave perturbations (SBgh = 0). The dis- 
tributions for <5-B s ;j/(52} = 50, and for 5B/Bq = 1.0 and 0.3 are 
presented with solid and dashed lines, respectively. The angular 
distributions for SB s f l / (B2) = 20, and SB/ Bo = 1.0 and 0.3 are 
shown with dash-dotted and dotted lines, respectively. 



to higher energies than that for 5B s h/(B2) = 50, which 
is caused by the large-scale component providing parti- 
cle scattering at higher energies (see, e.g., the spectrum 
for SBgh = in Fig. This is not the case for cither 
a flat spectrum (Fig. [Q) or a small-amplitude (Fig. EJ) 



of long- wave turbulence, since in superluminal shocks 
particle acceleration processes in the absence of largc- 
amplitud e long- wave magnetic field perturbations are in- 
efficient (|Begelman fc Kirkl Il990t iNiemiec &: Ostrowskil 
2004, 2006). 

Figure El presents particle angular distributions for 
superluminal shocks (V>i = 45°) with the Kolmogorov 
power spectrum of large-scale magnetic field perturba- 
tions. Note that the angular distributions for a large 
amplitude of the small-scale turbulence, 5B s h/ (B2) = 50, 
are closer to that for SB s h — (filled dots) than angu- 
lar distributions for the smaller amplitude, SB s h/ (B 2 ) = 
20. As the particle spectra, the distributions for 
o~B s h/ (B2) = 50 do not depend on the amplitude of the 
background field perturbations, nor on their wave-power 
spectrum (distributions for the flat power spectrum are 
not shown for presentation clarity). Similarity in shape 
of these distributions to the angular distribution derived 
without the 8B s h scattering term is caused by efficient 
particle scattering on the downstream short-wave pertur- 
bations at low particle energies, which enables a consid- 
erable fraction of them to return upstream of the shock. 
The difference in the distributions for the smaller short- 
wave amplitude, 8B s h/(B2) — 20, illustrates the effects 
the long- wave turbulence has on the particle acceleration 
process, as discussed above. 

A caveat is in order concerning the low-energy part 
of the simulated spectra. In our simulations with large- 
amplitude short-wave turbulence, the scattering ampli- 
tudes for low particle energies, E < 1 for the spectra 
with SB sh /(B 2 ) = 50 and E < 5 for SB sh /(B 2 ) = 80, 
violate the pitch-angle diffusion requirement A£7 <§; 1. 
This is because in our approach the time between subse- 
quent scatterings is fixed for a given particle energy, and 
the scattering angle scales linearly with 5B s h (see §2.2) 
providing large scattering amplitudes for large values of 
(B 2 ). Therefore, the simulated spectra may suffer a weak 
systematic flattening at low energies. 

3.2. Parallel Shocks 



As discussed bv lOstrowski fc Bednard l)2002|) . the par- 
ticle spectra with "universal" slope have been derived in 
conditions equivalent to a parallel ( , 0i = 0°) mean shock 
configuratio n. In our earlier work l|Niemiec fc Ostrowskil 
I2004L 120061) we have demonstrated that the spectra 
formed at parallel relativistic shocks depend substan- 
tially on background conditions, and features such as a 
relation between the spectral index and the turbulence 
amplitude or a cutoff formation in the case of ultrarel- 
ativistic shocks can be observed. In this section, we in- 
vestigate the role of shock-generated downstream turbu- 
lence in particle acceleration at parallel shocks. For that 
purpose, we have performed simulations for highly rel- 
ativistic shock waves with Lorentz factor 71 = 10 and 
30, as presented in Figures El and |H1 respectively. For 
parallel shocks the amplitude of the compressed back- 
ground field downstream of the shock, (B2}, is on aver- 
age much smaller than for the oblique shocks discussed 
in §3.1. Therefore, the numerical constraints arc largely 
relaxed, and we can use larger values of SB s h/ {B2), up 
to 300 or even 500. However, we would still slightly vi- 
olate the pitch-angle scattering approximation at E < 1 
for the spectra with 5B s h/ '(B 2 ) = 80 in Figs. Eb andEf), 
and 8B sh /(B 2 ) = 300 in Figs. Et andEJi. 
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Fig. 6. — Accelerated particle spectra at parallel ultrarelativistic shocks (71 = 10, = 0°). The parameters of the large-scale background 
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Figures^ and|HJ) show the results for shocks with 71 = 
10 and large-amplitude background field perturbations 
(SB /Bo = 1.0). The qualitative similarity of the spec- 
tral features observed in this case and for obliq ue superlu- 
minal shocks is apparent (see also lNiemiec &: Ostrowskil 



2006)— the spectra are convex and display cutoff's occur- 
ring much below E reSimax ~ 2ir/k m i n . This behavior 
proves that the large-amplitude long-wave perturbations 
locally form superluminal conditions at the shock, thus 
leading to a spectral cutoff as the pitch-angle scattering 
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term in the downstream shock-generated turbulence sub- 
stantially decreases with growing particle energy. Note, 
that in the case of a lower amplitude of the short-wave 
component, 8B s h/ (B2) = 20, the large-scale background 
turbulence with Kolmogorov distribution provides addi- 
tional particle scattering, allowing for the formation of 
a quasi-power-law spectrum in a slightly wider energy 
range than for the larger 8B s h/{B2) (Fig. HJ)), again in 
correspondence to the analogous result for the superlu- 
minal shocks (Fig. 

For the smaller amplitude of the large-scale magnetic 
field perturbations, 5B/B = 0.3, the role of the long- 
wave perturbations is less significant and the scattering 
on the short-wave turbulence can dominate up to the en- 
ergies higher than the upper resonance energy E res>max . 
Though being gradually weaker with growing particle 
energy, in this shown in Figure HJ:, pitch-angle 

diffusion can dominate over scattering by the long-wave 
perturbations without a limit for particle energy below 
E re s,max- The scattering mean free paths along the mean 
background magnetic field (shock normal) and, corre- 
spondingly, the mean acceleration timescale also increase 
with particle energy, but the particle spectra retain a 
power-law form up to the highest energies E w E res ^ max 
studied in our simulations. However, the spectra are 
steeper than the expected "universal" spectrum, a > a u , 
and only weakly dependent on 8B s h ■ To understand this 
feature, we have performed additional test simulations 
with the short-wave component introduced both down- 
stream and upstream of the shock (with constant up- 
stream scattering amplitude AQ(E) up — const <C 1/71)- 
This case, equivalent to the pure pitch-angle diffusion 
model, should lead to the formation of the "universal" 
spectrum and, indeed, such the spectrum is produced 
in our simulations (described as "scatt up" in Fig. &p). 
Thus, the slightly steeper spectra obtained in the mod- 
els with only a downstream short-wave component can be 
understood solely as the result of the upstream long- wave 
turbulence providing a different distribution of pitch- 
angle perturbations to the particle orbits than the purely 
diffusive process (see Fig. 0. This allows for the gener- 
ation of a power-law particle spectrum, but with some- 
what steeper index than the "universal" value. 

With the Kolmogorov background wave spectrum (Fig. 
EK) there is enough power in long waves to influence the 
spectrum at high energies, leading to the noticable steep- 
ening of the particle distributions below E res ^ max . How- 
ever, also in this case the model with the short-wave field 
perturbations imposed upstream of the shock yields the 
"universal" power-law spectrum. In Figure we com- 
pare the angular distributions of particles at the shock 
in the cases of only downstream small-scale turbulence 
and with pitch-angle diffusion also upstream of the shock. 

The spectra for weakly perturbed background mag- 
netic fields, 8B/B n = 0.1, presented in Figures Els and 
EK show a new characteristic feature, namely a spec- 
tral bump besides the power-law spectrum at lower 
particle energies. The bump shape closely resembles 
the distribution of particles accelerated without the 
shock-generated field comp onent (see the discussion in 
iNiemiec fc Ostrows ki 2006), and its presence indicates 
that at these energies the short-wave perturbations be- 
come negligible compared to the background long-wave 
component. Note, that 8B s h is scaled with respect to 



: <5B sh /<B 2 >=100 


- 


; 6B/B =0.3 j\ 




- 7i = 10 / i 


- 


— ~~t— / \ 


- 


F(k)~k- 5/3 : '7\ \\ 


- 


! / y \ 


\ *\ c/^/tt ft 1 in 

\ \ 




\ ^ - 


1:1 


\ t 


■>/ 


• V- 




■ \ v ' ~ 

\ \ v. 


, 1 , , , , 1 


\. V. 



-1.0 -0.5 0.0 0.5 1.0 

cos(e) 



Fig. 7. — Particle angular distributions at parallel shocks (71 = 
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(solid and dashed line) or the Kolmogorov (dash-dotted and dotted 
line) wave power spectrum of magnetic field perturbations. The 
distributions indicated as " scatt up" are derived in a model with 
short-wave turbulence imposed both downstream and upstream of 
the shock. 



(i? 2 ), which is relatively small in comparison to the cases 
with larger 8B/B discussed above. 

We have also performed a series of simulations for 
shocks with a larger Lorentz factor of 71 = 30. The par- 
ticle spectra presented in Figure El for 8B s h/{B2) = 100 
show a similarity of spectral features to those observed 
for shocks with 71 = 10. However, there is an approx- 
imate scaling between the two cases, with characteris- 
tic spectral features appearing at respectively lower am- 
plitudes of the long-wave magnetic field perturbations 
8B/Bq for shocks with 71 = 30 than in the spectra for 
ji = 10. Such a scaling may result from the stronger 
compression of these perturbations, as measured between 
the two plasma rest frames, and also from the larger par- 
ticle anisotropy involved at the higher-7 shock. 

4. CONCLUDING REMARKS 

The present paper conclu des our studies 
ijNiemiec fc Qstrowskil 12004 I2006D of the first-order 
Fermi acceleration mechanism acting at relativistic and 
ultrarelativistic shock waves. We have attempted to 
consider as realistic models as possible for the perturbed 
magnetic field structures at the shock, which allow us to 
study all the field characteristics important for particle 
acceleration. In particular, we have investigated the 
dependence of particle spectra on the mean magnetic 
field inclination with respect to the shock normal and 
on the power spectra of magnetic field perturbations, 
including the long-wave background turbulence and the 
short-wave turbulence generated at the shock. 

In the present p a per w e append the model of 
INiemiec fc Ostrows ki (2006J) with downstream large- 
amplitude small-scale MHD turbulence component, anal- 
ogous to those seen in PIC simulations of collisionless 
relativistic shocks. The form of these additional pertur- 
bations is arbitrarily chosen to consist of a very short 
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sinusoidal waves, that form (static) isotropic highly non- 
linear turbulence. We have considered a wide variety of 
ultrarelativistic shock configurations for which, as pre- 
sented in §3, a number of different spectral features can 
be observed. We show that with growing particle en- 
ergy the role of short-wave (A s /j <C r g (E)) magnetic 
field perturbations decreases, and spectra at high ener- 
gies are shaped only by the mean magnetic field and the 
long- wave perturbations. Thus, in oblique superlumi- 
nal shocks concave particle spectra with cutoffs can be 
generated even for 6B s h ^> (B2). Extended power-law 
particle distributions can be formed in parallel shocks 
propagating in a medium with low-amplitude long-wave 
field perturbations, but the spectral indices obtained are 
larger than the "universal" spectral index a u m 4.2, 
that is widely considered in the literature. The only 
case in which we have been able to obtain spectra with 
a = a u involved the unphysical assumption, that large- 
amplitude short-wave field perturbations also exist up- 
stream of the shock. 

In summary, the accelerated particle spectral distribu- 
tions obtained in this work and in our previous studies 
l|Niemiec fc Ostrowski 12004 l200fA generally differ from 
the spectra of relativistic electrons (power-laws with the 
spectral indices close to a u ) inferred from modeling the 
electromagnetic emission spectra of astrophysical sources 
hosting relativistic shocks (e.g., hot spots in extragalactic 
radio sources, quasar jets, gamma-ray burst afterglows). 
Our results thus provide a strong argument against con- 
sidering the first-order Fermi shock acceleration as the 
main mechanism producing the observed radiating elec- 



trons. In our opinion, other processes must be invoked 
to explain the observed emission spectra, for example 
second-order Fermi processes acting in the regions of rela- 
tivistic MHD turbulence downstream of the shock, which 
have hardly e ver been discussed for r elativistic condi- 
tions till now l|Virtanen fc Vairlioll2005l) . or the collision- 
less plasma processes that have been studied in numerous 
PIC simulations (e.g..|Hoshino et al.ll992tlHededal et alJ 
120041 iNishikawa et alJl2005bl) . or other no n-standard ac- 
celeration proc esses like those discuss ed bv lSternl ll2003fl 
for electrons or lDerishev et all l)2003fl for both electrons 
and protons. 

Our simulations also show, that relativistic shocks, 
being essentially always superluminal, possibly gener- 
ate accelerated particle distributions with cutoffs below 
either the maximum resonance energy enabled by the 
high- amplitude background turbulence (r g (E cuto ff) < 
A(£J res>max )), or approximately at the energy of the 
compressed background plasma ions E cuto j f ~ 7imjC 2 
(where m, is a mass of the heaviest ions present in the 
background medium). Thus, in conclusion, we main- 
tain our opinion from th e previous publications (see also 
iBegelman fc Kirklll99fl) that relativistic shocks are not 
promising sites as possible sources of ultra-high-energy 
nuclei registered by the air shower experiments. 
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